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Jet velocity in SS433: its anti-correlation with precession-cone 
angle and dependence on orbital phase 

Katherine M. Blundell 1 and Michael G. Bowler 1 
ABSTRACT 

We present a re-analysis of the optical spectroscopic data on SS 433 from the 
last quarter-century and demonstrate that these data alone contain systematic 
and identifiable deviations from the traditional kinematic model for the jets: vari- 



>- ■ ations in speed, which agree with our analysis of recent radio data; in precession- 



cone angle and in phase. We present a simple technique for separating out the 



in 

jet speed from the angular properties of the jet axis, assuming only that the 
jets are symmetric. With this technique, the archival optical data reveal that 
the variations in jet speed and in precession-cone angle are anti-correlated in the 
sense that when faster jet bolides are ejected the cone opening angle is smaller. 
CLi 1 We also find speed oscillations as a function of orbital phase. 
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Subject headings: stars: binaries: stars: individual: (SS 433) 



^ ' 1. Introduction 



In a recent paper (Blundell & Bowler 2004) we presented the deepest yet radio image 
of SS433, which revealed an historical record over two complete precession periods of the 
geometry of the jets. Detailed analysis of this image revealed systematic deviations from the 
standard kinematic model (Margon 1984; Eikenberry et al 2001). Variations in jet speed, 
lasting for as long as tens of days, were needed to match the detailed structure of each jet. 
Remarkably, these variations in speed were equal, matching the two jets simultaneously. 

The Doppler residuals to the kinematic model show little variation with the precessional 
phase of the jets and this observation rules out variations in jet speed alone as the source of 
the residuals (e.g. Katz & Piran 1982; Eikenberry et al 2001). Very little phase variation 
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is obtained if the pointing angle jitters (Katz & Piran 1982) but there is no evidence ex- 
cluding symmetric speed variations of the magnitude reported in Blundell & Bowler (2004) 
superposed on pointing jitter, Figl. Thus our findings from the radio image led us to re- 
analyse the archival optical data. G. Collins II and S. Eikenberry (with kind permission of 
B. Margon) made available their compiled datasets, published in Collins & Scher (2000) and 
used in Eikenberry et al (2001). We use the Collins' compilation (available at http://www- 
astro.physics.ox.ac.uk/~kmb/ss433/) because of its higher quoted precision, but very similar 
results are obtained from Margon's. 



2. Speed and angular variations from the optical data 

If the precessing jet axis of SS 433 traces out a cone of semi-angle 9 about a line which is 
oriented at an angle i to our line-of-sight with jet velocity (3 in units of c (7 = (1 — /3 2 ) -1 / 2 ), 
the redshifts measured from the west jet (z + ) and the east jet are given, if the jets are 
symmetric, by: 

z± = —1 +7[1 ± (3 sin 9 sin i cos ± (3 cos 9 cos i] , (1) 

where is the phase of the precession cycle (see http:/ /www-astro. physics. ox.ac.uk/~kmb/ss433/). 
Addition of z + and z^ in Eqn 1 gives an expression relating the observed redshifts to the jet 
speed independently of any angular variation. Re-arrangement gives 

1/2 

(2) 

The quantity z + + z_ fluctuates very substantially (see Fig 2). Subtraction of the expressions 
for z + and Z- gives the angular properties a of the orientation of the jet axis, with the speed 
divided out using Eqn 2: 

a = + = sin 6* shu cos0 + cos 9 cost, (3) 

Fluctuations in (3, 9 and (b are predominantly symmetric, as described by Eqns 1-3, when 
fluctuations in z + + z_ represent symmetric fluctuations in speed. (The velocity variations 
in the two radio jets (Blundell & Bowler 2004) are highly symmetric; the standard deviation 
on the difference in the speeds is less than 0.004 c and on the common velocity 0.014 c.) 

The disadvantages of the variables s = z + + Z- and a (Eqn 3) are that their interpreta- 
tion is simple only for perfect symmetry and that they may only be used for the 395 out of 
486 observations which record a simultaneous pair. Their merits are exemplified by Fourier 



P 



1 - 



1 + 



z + + z_ 



1 -2 



- 3- 



analyses of the time distributions. We used the algorithm of Roberts et al (1987) which 
accounts for the uneven time-sampling of the data. The angular data a clearly revealed pe- 
riodicities corresponding to the nodding of the precession axis (Katz et al 1982; Newsom & 
Collins 1982; Collins & Scher 2002) and the 162-day precession period, clearly seen in Fig 3a. 
There is no periodicity in the speed data s (Fig 3b) common to the angular data, consistent 
with perfect symmetry. The speed data also indicate a periodicity at 13.08 days (the period- 
icity at 12.58 days matches a beat with Earth's orbital period 365 days); to investigate this, 
we folded the data over 13.08 days in 20 phase bins, and in each bin the mean speed ((3, from 
Eqn2) was derived. Fig 4 shows a clear sinusoidal oscillation with orbital phase. The rms 
variation in speed which oscillates with orbital phase is smaller by a factor of three than the 
overall speed dispersion. This oscillation with amplitude 2000 km s -1 may be because the 
speed with which the bolides are ejected is a function of orbital phase, but the excursions in 
Fig 4 could also be interpreted as due to orbital motion; in that case SS433's orbital speed 
is ~ 400kms~\ 



We fitted Collins' dataset with the kinematic model, including nodding. From our fit, 
we derived model redshift pairs and hence the variables s and a. Subtraction of these model 
variables from those constructed from the data gave residuals As in s and Aa in a. The f3 
variation is shown in Fig 2c. The standard deviation of this histogram is 0.013, in excellent 
agreement with the result from our radio image (0.014). Examples of the residuals in s and 
in the angular variable a are plotted in Fig 2; the variations in speed and angular residuals 
are anti-correlated. From Eqnsl-3, maintaining the assumption of symmetry: 



3. Anti-correlated deviations in jet speed and 9 



As 2 = 4/3 V A/3 2 , 



(4) 



Aa 2 = (cos 6* sin i cos — sin 6> cosi) 2 A# 2 + (sin6 l sinisin0) 2 A0 2 
— 2sin#sinisin0(cos#shucos</> — sin 9 cosi) A6A0, 



(5) 



Aa As 



2/37 3 [(cos6 ) sini cos0 — sin 9 cos i) Af3 A9 
sin 9 sin i sin 0A/3A0] , 



(6) 



where A/3, A9 and A0 represent the variations in (3, 9 and respectively. Averaging over 
many cycles for any given value of the phase yields the averages (Aa 2 ), (As 2 ), (AaAs) 
as functions of 0, in terms of the parameters (A/3 2 ), (A9 2 ), (A/3 A9) and so on. The fit 
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to (A(3 A9) is shown in Fig 5. The parameters are given in Table 1; the global y^lNDF is 
34.8/24, where NDF is the number of degrees of freedom. 

Fig 5 shows that the quantity (AaAs) has an almost pure cosinusoidal variation with 
phase, as given by the first term on the right hand side of Eqn6. This shape is the unique 
signature of a correlation between variations in f3 and 9. (Af3 2 ) shows no correlation with 
and (Aa 2 ) very little. The latter requires fluctuations in both 9 and in 0. We remark 
that (A(z + — Z-) As) does not show a strong correlation with 0; nor should it, using the 
parameters from Table 1. Removing the varying speed from z + — z_ was crucial in revealing 
this correlation in (AaAs). 

4. The redshift residual plot 

Consider the plane of redshift residuals, as in figure 5 of Eikenberry et al (2001), with 
symmetric excursions from the kinematic model in /3, 9 and 0. Comparison of their figure 5 
(similar to our Fig6f) with our Figs 6a and 6b requires the presence of both angular variations 
(to spread the points along the line y = —x) and velocity variations (to spread the points 
perpendicular to this line — note that even their quoted redshift measurement error of 0.003, 
likely an over-estimate, will not account for this breadth). Inclusion of all these variations, 
correlated as in Table 1, gives Fig6d which resembles that from the data (Fig6f). The 
simulations in Fig6d take no account of the (stochastic) duration of the excursions; the 
duration of variations in Fig 6e were drawn from a gaussian with half-width 2 days. 

Thus Figs 1 and 6 establish the consistency of the optical data with perfect symmetry 
and with speed fluctuations whose magnitude is in excellent agreement with those found in 
the radio image. In addition, the slope of Fig6e is —0.765 ± 0.034, and that for Fig6f from 
the Collins data set is —0.786 and not —1; this long standing curiosity is explained by the 
physics from the § 3 fit. 

5. Other assumptions 

If the jet speed were constant, the residuals to z + and z_ would, for the case of strictly 
antiparallel jets, be correlated as in Fig 6a; to spread the distribution of points perpendicular 
to this line it is necessary to allow some independence in the pointing of the two jets. The 
variation with precessional phase of the quantities (As 2 ), (Aa 2 ) and (AaAs) is then almost 
as well described (x 2 /NDF = 36.3/24) as by our fit of §3. Such a model would not be able 
to explain the radio jet morphology (Blundell & Bowler 2004) and this fit required angular 
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fluctuations breaking symmetry to have rms values ~ 1/2 of those preserving symmetry. 
Angular jitter alone cannot account for the slope of Az_ versus Az + (Fig 6) differing from 
— 1, unless the jitter in the East jet is systematically smaller than in the West jet. Velocity 
variation breaks symmetry in the Doppler shifts and accounts naturally for this observation. 
A better fit than either (x 2 /NDF = 25.5/21) was achieved by allowing some symmetry 
breaking angular fluctuations in addition to symmetric velocity fluctuations: in this case the 
rms symmetry breaking fluctuations were ~ 1/4 of those preserving symmetry. In all cases 
the rms fluctuations in 6 were ~ 1/3 of the rms fluctuations in 0, as would be expected for 
pointing angle fluctuations described by Katz & Piran (1982) as "isotropic". 



6. Concluding remarks 

Archival optical spectroscopic data on SS433 reveal variations in jet speed, in cone 
opening angle, and in the phase of the precession. These appear in the plane of redshift 
residuals (Fig 6) and through the new (symmetry dependent) technique of combining simul- 
taneously observed redshift pairs for the speed-only (s) and angular-only (a) characteristics 
(Fig 2). The velocity variations ~ 0.014 c are strongly anticorrelated with cone angle in the 
sense that when faster bolides are ejected the cone angle is smaller. We also found smaller 
amplitude sinusoidal oscillations in speed as a function of orbital phase. If this is due to 
ejection speed, perhaps the orbit of the binary is eccentric. If these 13.08-day oscillations are 
orbital Doppler shifts the orbital velocity is ~ 400 km s" 1 (twice that inferred by Crampton 
& Hutchings (1981); Fabrika & Bychkova (1990)). If this were the case, the mass of the 
companion to SS433 would be > 86 M and if the mass fraction were 0.1, then the mass of 
the companion would be > 100 M Q . Such masses would be hardly consistent with an A-type 
companion (Gies et al 2002; Charles et al 2004). 

K.M.B. thanks the Royal Society for a University Research Fellowship. It is a pleasure 
to thank Avinash Deshpande, James Binney & Philipp Podsiadlowski for helpful discussions. 
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Table 1. Fits to excursions from the standard kinematic model (incorporating nodding). 

(3 is in units of c, and 9 and (f> are in radians. 



Quantity 


fitted value 


quantity 


derived value 


(A/? 2 ) 


1.67 ±0.18 x lCT 4 


rms speed variation 


0.0129c 


(A9 2 ) 


2.24 ±0.35 x 10~ 3 


rms 9 variation 


2.71 deg 


(A0 2 > 


1.31 ±0.26 x 10~ 2 


rms variation 


2.96 days 


{A/3 AO) 


-3.81 ±0.52 x 10~ 4 






(A/3A0) 


1.70 ±0.80 x 10- 4 






(A#A0) 


indistinguishable from zero 
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Fig. 1. — Demonstration that the dispersion in redshift with phase is unaffected by velocity 
variations as long as they are accompanied by pointing variations: (a) random variations in 
velocity drawn from a gaussian of half-width 0.013c, (b) random variations in 9 and in phase 
drawn from gaussian distributions of 2.7 deg and 2.95 days respectively and (c) uncorrelated 
variations in all three of the quantities above. 



- 9- 




-0.03 1 ' 1 ' ' ' ' 1 ' 1 ' ' 1 

4000 4010 4020 4030 4040 4050 4060 

Modified Julian Day 



Fig. 2. — (a) and (b) examples of angular and speed residuals versus time. The variations 
suggest anti-correlation, §3. (c) Distribution of speeds from the data compiled by Collins. 
The mean is indicated by the vertical red line, and the standard deviation in (3 is 0.013. 
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Fig. 3. — Fourier transform of data as described in §2. (a) The differences of all redshift 
pairs, divided by 2/?7 (see Eqn 3) which, if the jets are symmetric, depend only on the angular 
properties of the jet. (b) The sums of all redshift pairs which, if the jets are symmetric, 
only depend on the speed of the jet. 
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Fig. 4. — Speed data (Eqn2) folded over the orbital period of 13.08 days showing a clear 
sinusoid, mean 0.2581 ± 0.0005, amplitude 0.0066 ± 0.0007, phase offset with respect to 
optical ephemeris 2.17 ± 0.11 rad. 
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Fig. 5. — The mean value of (Aa x As) in 162-day period phase bins, showing the best fit 
to Eqn6 averaged. There are two free parameters in the fit shown as a solid line, (A/3 A9) 
and (A/3 A0) . The dotted line shows the fit with fixed velocity but angular asymmetry. 
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Fig. 6. — The effect on redshift residuals in each jet for deviations from the kinematic model 
in 9, (3 and 0, for each date on which there is a redshift pair (replicating the sampling 
function), (a) Variations only in 9, drawn from a gaussian distribution with half-width 
2.7 degrees. The form of this plot, points lying only on the y = —x line, is identical for 
variations only in 0, or for variations in both and 9. (b) Speed-only variations (on a 
given day (3 is drawn from a gaussian of half-width 0.013 — Fig 2c). (c) Variations in j3 
anti- correlated with those in 9: on a given synthesized observation date the same randomly- 
drawn number from a gaussian is scaled by 0.013 for the speed variation and by —2.7 degrees 
for the 9 variation, (d) As (c), but with 9 partially anti-correlated with (3 and uncorrelated 
fluctuations (as Table 1). (e) As (d), but the duration of the /3-variations, and that of 
the correlated ^-component, is drawn from a gaussian of half-width 2 days, (f) Observed 
residuals. 



